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Abstract 



The edge state theory of a class of symmetric double-layer quantum Hall 
systems with interlayer electron tunneling reduces to the sum of a free field 
theory and a field theory of a chiral Bose field with a self-interaction of the 
sine-Gordon form. We argue that the perturbative renormalization group flow 
of this chiral sine-Gordon theory is distinct from the standard (non-chiral) 
sine-Gordon theory, contrary to a previous assertion by Renn, and that the 
theory is manifestly sensible only at a discrete set of values of the inverse 
period of the cosine interaction (/?). We obtain exact solutions for the spectra 
and correlation functions of the chiral sine-Gordon theory at the two values 
of (3 at which electron tunneling in bilayers is not irrelevant. Of these, the 
marginal case (p 1 = 4) is of greatest interest: the spectrum of the interacting 
theory is that of two Majorana fermions with different, dynamically generated, 
velocities. For the experimentally observed bilayer 331 state at filling factor 
1/2, this implies the trifurcation of electrons added to the edge. We also 
present a method for fermionizing the theory at the discrete points (/3 2 G Z + ) 
by the introduction of auxiliary degrees of freedom that could prove useful in 
other problems involving quantum Hall multi-layers. 
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I. INTRODUCTION 



This paper is concerned with a chiral version of the celebrated sine-Gordon (SG) field 
theory in one dimension, in which the standard kinetic term is replaced by the kinetic term 
for a chiral Bose field. More precisely, we are concerned with the Hamiltonian for a periodic 
boson of radius R — 1/(3, 



7-L x sg14>} = I dx 



(1) 



An 

supplemented by the equal-time commutation relation, 

[4>(x), 4>(x')] = in sgn(x — x'), (2) 

appropriate for a chiral Bose field. The Hamiltonian is suggestive of the standard sine- 
Gordon theory, but the resemblance is misleading and the physics, as we shall see, is very 
different on account of the chiral nature of the field. Indeed, a previous analysis of this 
problem by Renn [|IJ is in error precisely because of his neglect of the chiral constraint. We 
will comment further on the difference between the two theories at several points in the 
main text. 

The chief interest of the chiral sine-Gordon (%SG) theory is that it arises naturally, as we 
detail below, in the edge theory of double-layer quantum Hall systems when it is necessary 
to take account of weak interlayer-tunneling within Wen's chiral bosonic description of the 
edge dynamics. As such it is the simplest of a family of interacting chiral theories in one 
dimension that describe the dynamics of multicomponent systems. (Such systems yield 
non-chiral problems as well, in cases where there are counter-propagating modes, and they 
have attracted considerable attention J2|-|J.) These theories are much more constrained 
than generic non-chiral theories, but they are by no means trivial and they exhibit unusual 
connections between renormalization group flows and spectral rearrangements. In a sense, 
they are intermediate in complexity between quantum impurity problems, which can be 
mapped to chiral problems with interactions localized at a point, and generic one dimensional 
field theories, and seem worth studying on this account as well. Dijkgraaf has explored some 
general properties of chiral theories from a conformal field theory viewpoint ||, and several 
authors have considered particular examples such as the chiral Potts model || and the 
chiral Schwinger model [0]. Finally, the study of chiral field theories is not without interest 
in another condensed matter context: Ho and Coleman have recently motivated the study 
of interacting chiral Majorana theories by appealing to the phenomenology of the cuprates 
and have presented solutions of several models by an interesting "fermionic bootstrap" ||. 

We begin with a review of the bosonic formulation of the edge theory of double-layer 
quantum Hall systems with interlayer tunneling (Section |II|). For systems with one mode 
per layer we show that for clean, symmetric realizations, the edge theory separates into the 
sum of a free boson Hamiltonian and the %SG Hamiltonian with (3 2 G Z^, en (the "bosonic" 
sequence). We argue that the x^G theory is manifestly sensible only for integer values of 
p 1 . Layer asymmetry and disorder, which lead to a more complicated dynamics, will be 
discussed elsewhere ||. Next we consider the renormalization group flow of the %SG theory, 
perturbatively in A, and argue that it differs from the usual SG theory (Section |nj). We 
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argue that (3 is not renormalized and for p 2 = 4 there is a line of fixed points parameterized 
by A. In Section [TV] we obtain exact solutions for the spectrum and correlation functions of 
the xSG theory for p 2 = 2 and p 2 = 4 which are the two members of the bosonic sequence 
at which the interaction is not irrelevant. At (3 2 = 2 we show by "re-bosonization" , that the 
spectrum of the interacting theory is, in the thermodynamic limit, still that of a free chiral 
boson. Perhaps of greatest interest is our solution by fermionization at p 2 = 4, relevant to 
the experimentally observed 331 state, where the xSG theory is shown to be equivalent to 
a theory with two chiral Majorana fermions with dynamically generated distinct velocities. 
Next we present a method for fermionizing the %SG theory at p 2 G Z + by the introduction 
of auxiliary bosons (Section |V|). We show that for p 2 = 2 and p 2 = 4 the fermionized theory 
with auxiliary degrees of freedom can be solved exactly, and after projection onto the %SG 
Hilbert space we recover the results of Section |IV|. In Section [VI] we briefly consider the 
%SG theory at the points f3 2 G Z^ dd (the "fermionic" sequence), which do not correspond 
to any double-layer system. We comment briefly on the difficulties in understanding our 
results from the kinds of semiclassical considerations that are so useful in the standard sine- 
Gordon problem in Section VII. We conclude with a discussion (Section VIII ) and several 
Appendices giving detailed calculations in support of statements made in the text. 



II. EDGE THEORY OF DOUBLE-LAYER SYSTEMS 

A. Double-layer Systems without Tunneling 

We are interested in a system consisting of two parallel layers of 2DEGs (two-dimensional 
electron gases) in a strong perpendicular magnetic field, with a confining potential that 
restricts the electrons in each layer to a region with the topology of a disc. The geometry 
is sketched in Fig. p]. The simplest Abelian quantum Hall states in double-layer systems 
are well described by the two-component generalization of the Laughlin wavefunction first 
introduced by Halperin fli~0|| : 



*W,n({*ta}) = II - zip) m (z2a ~ z 2 p) m ' ]J{z la - z w ) n e~ ^ |2/4 , (3) 

a<(3 a, p 

where z ia is the complex coordinate of electron a in layer i, and we work in the symmetric 
gauge. Here m and ml are positive odd integers characterizing the intralayer correlations 
and n is a non-negative integer specifying interlayer correlations. The filling factor in each 



layer is 11 



u\\ 1 / m' — n i 

(4) 



i>2 ) (mm' — n 2 ) \ m — n 
Following convention we refer to these states by the shorthand mm'n. 
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FIGURES 




FIG. 1. The overall geometry of a double-layer quantum Hall system, with edge modes in 
both layers (ui) propagating in the same direction. The chiral sine-Gordon theory describes the 
dynamics of the neutral mode in the presence of uniform interedge electron tunneling. 



The gapless excitations of the system are confined to the edges of the droplets, and for 
unreconstructed edges they are described at low-energies by the effective Hamiltonian |L^] 

rL/2 1 

Ho = / dx—Vij :d x Uid x uf., (5) 

J — L J '2 4rVT 

where we have assumed the long-range part of the Coulomb interaction is screened. Here x 
is the coordinate along the edge of length L and normal ordering will be denned below. V 
is a symmetric matrix whose elements depend on the details of the confining potentials and 
the interactions at the edge. We require V to be positive definite so that the Hamiltonian 
is bounded from below. We can parameterize this matrix by 

V — (j ^ J , where g 2 < vv'. (6) 

The Ui(t,x) are chiral Bose fields with compactification radius R = 1 (i.e., u ~ u + 2ir), 
which satisfy the equal-time commutation relations: 

[ui(t, x),Uj(t, x')} = inKij sgn(x — x'), (7) 

where K is a symmetric, integer-valued matrix which characterizes the topological properties 
of the system and is completely determined by the exponents in the bulk wavefunction [13| 



n m 



K=C !,)■ (8) 



From here on, any expression in which the time arguments of the operators are suppressed 
is to be understood to hold at equal times. Eq. (|7|) is valid in the limit of large L; if we 
retained the full L dependence the r.h.s. would be periodic (mod 2iriKij) in (x — x') with 
period L. 

Since the chiral Bose fields are compact, we can introduce the integer- valued Hermitian 
topological charge operators 

1 r L / 2 

M = — y ^^dxd x Ui(x), (9) 
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which by Eq. (0) obey: 



[Afi, Uj (x)} = iKij, [M, A/}] = 0. 



(10) 



For a chiral boson with a radius R ^ 1, the coefficient of the integral in Eq. (|9|) would be 
1/(2ttR). The electric charge density on the edge of layer i, Pi{x), and the total charge, Qi, 
are related to the corresponding topological quantities by: 



(11) 



The vertex operators which create an edge quasiparticle or an edge electron in layer i are, 
respectively, 



qp,i\ 



*Ux) oc e~ lUl{x) . 



(12) 



These operators change the topological and electrical charges according to 



l 3 qp,J ' 

K -vpt 



t 



(13) 



The physical meaning of the edge quasiparticle creation operator is as follows. When a 
quasihole is created in the bulk fluid {e.g., by the insertion of a flux quantum), the charge 
depleted from the vicinity of the quasihole accumulates at the boundary of the droplet. In the 
somewhat artificial theoretical construction in which the "bulk" and "edge" of the droplet 
are considered distinct, this accumulation appears as the addition of a fractional amount of 
charge to the edge theory, which is accomplished by the edge quasiparticle creation operator. 



If we consider the action of the vertex operators ([12]) in a two-dimensional integer lattice 
in the A/iA/2-plane, the quasiparticle creation operators can be identified with the basis 



vectors of this lattice: \l/ 



t 



electron creation operators, ^\ 



[m, n 



qp,2 



(0, 1). The vectors corresponding to the 
(n, m'), are the basis vectors of an "electron 



sublattice" consisting of all states connected to one another by adding or removing only 
electrons from the edges 0]. The number of distinct electron sublattices is equal to the 
volume of the unit cell spanned by the electron vectors: mm' — n 2 . We shall refer to these as 
different sectors of the edge theory. Note that the discrete spectrum of the topological charges 
also places restrictions on the allowed vertex operators, as we describe below. An example of 
the topological charge lattice is shown in Fig. || for the 331 state {i.e., m — m! — 3,n = 1). 
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FIG. 2. The topological charge lattice in the A/iA/2-plane for the double-layer 331 state. The 
crosses denote points in the electron sublattice (i.e., sector) containing the origin. The unit cell of 
this sublattice is enclosed by the dashed line. 



From Eqns. fllPf) and (|T^) we find 



— 2iriAfi —i-KKi, 



+ L) = *l p ,(x)e- 2 ™ K ^e-^\ (14) 

Since A/j, G Z the boundary conditions on the electron operators are independent of the 
values of the topological charges. Within the sector containing the origin, Mi = 0, the same 
statement holds for the quasiparticle operators. However, in other sectors (e.g., the state 
with a single quasiparticle in layer 1, Af\ = 1, A/2 = 0) the quasiparticle operator acquires an 
additional phase as it is taken around the edge. Since creating a quasiparticle at the edge 
requires the creation of a quasihole in the bulk as noted above, this additional phase can be 
interpreted as the Aharonov-Bohm phase due to the flux creating the bulk quasihole |[L5|| . 

We are interested in the case where the two edge modes have the same chirality. Since 
V is positive definite, the chirality of the edge bosons is determined by the signs of the 
eigenvalues of K. Therefore we shall restrict our discussion to the "maximally chiral" case 
in which both eigenvalues of K are positive, which requires n 2 < mm'. (The requirement 
that the filling factor in each layer be positive then gives n < m, m'.) In this case, the scaling 
dimensions of vertex operators are independent of V ||16|1 . Furthermore we shall consider 



only the symmetric case: m = m' and v = v'. Our result in Section [QJ for the lowest-order 
perturbative RG flow of the tunneling amplitude applies for the case v ^ v' and can easily 
be extended to states with m 7^ m', but the mapping to the chiral sine-Gordon theory exists 
only in the symmetric case, as we demonstrate below. 

The Hamiltonian, Eq. (|5|), can be simplified by performing an orthogonal transformation 
on the chiral Bose fields to diagonalize K and V, followed by a rescaling to bring K to the 
identity. The combined transformation to the new fields C , <p n is: 



ui \ 1 / \Jm + n —\/m — n 
U2J y/2 V V m + n V m - n 



(15) 
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We will refer to <p c , n as the charged and neutral edge modes, respectively, since the topo- 
logical charge of C is proportional to the total electric charge on the two edges while the 
topological charge of 4> n is proportional to the electric charge difference between the edges. 
When we write i; the index i runs over the two values i = c, n. In terms of the new fields, 
the commutators are now independent and conventionally normalized: 

[(f>i{x), (f>i(x')] = iirdij sgn(x - x). (16) 

The mode expansions of the rotated chiral bosons are 

9 -i oo i 

<Mx) = — N&x - - X i + E ^e~ Kal \e- %Kx h\ r + e ik ^b ir ), (17) 

where k r = 2nr/L, R { are the compactification radii of the rotated fields, iV; are the new 
topological charges, Xi are the conjugate Hermitian phase operators, b ir are bosonic anni- 
hilation operators, and a is a short-distance cutoff. The only nonvanishing commutators 
are 

[b ir , bj s ] = 5 ij( 5 rs , [xuN^iSiy (18) 

The normal ordering introduced in Eq. ([|) is defined with respect to the bosonic oscillators 
bi r ,b\ r . Up to this point we have neglected normal ordering when writing vertex operators. 
Using Eq. (|T7|) one can show that normal ordering produces a multiplicative renormalization 
of vertex operators which to lowest order in a/L gives 

\2na) ' V ' 

where A(a) = a 2 /2 is the scaling dimension of the operator, and also its conformal spin 
since 0i is a chiral field. 

Using Eqns. flT5| ) and ( |TTD we can express the Hamiltonian (^|) in terms of the modes of 
the uncoupled fields: 

r L / 2 r 1 1 ttv- 2ir 00 

H = <l.r —r : :(<),,:>, f: = -i (R^) 2 + — £ w^An (20) 

J-L/2 L47T J L L ~^ 



where the velocities of the eigenmodes are v c>n = (m±n)(v±g). The structure of the Hilbert 
space is now clear. There are an infinite number of states \N c ,N n ) which are eigenstates of 
the topological charge operators and which are annihilated by b ir for i = c, n, and r > 0. This 
infinite set of oscillator vacua is divided into a finite number of sectors as defined above. The 
(zero-energy) ground state is |0, 0). Note that because Eq. ( |l"5l) defines a homogeneous linear 
transformation, the states |A^ c ,iV n ) are also vacua for the oscillator modes of the original 
bosons (ui). 

The values of the radii R\ are determined by the requirement that the new winding 
number operators N- Y have integer eigenvalues. From the zero-mode piece of the inverse of 
the transformation fll5T) we find 
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R C N C = 1 (M+AT 2 ), R n N Q = 1 (7V 2 -M)- (21) 
'2(m + n) \ 2{m — n) 



Since A/i G Z we see that the choice _R C = 1/ ^2{m + n), _R n = 1/ \j2(m — n) restricts 
N c — J\fx + A/2, A^ n = A2 — A/i to be integers. From this relation we also see that iV c and 
A" n are not independent since they must satisfy the "gluing condition" N c + N n E Z cven in 
order to describe the original A/1A/2 integer lattice. If one considers the electron sublattice 
containing the ground state (i.e., the origin A/i = N{ — 0), then one may choose a unit cell 
such that the m? — n 2 points in distinct electron sublattices are labeled by 

{{N c , N n )\N n = 0,-l,...-(m-n-l);N c - F Nn = 0, 2, . . . 2(m + n - 1)}, (22) 

where the offset _F/v n = [1 — (— l) Nn ]/2 is zero if N n is even and one if N n is odd. From this 
it is clear that the m 2 — n 2 distinct sectors of the full edge theory correspond to f n = m — n 
sectors for the neutral boson and f c = m + n sectors for the charged boson. The requirement 
that the vertex operator e~ %OL ^ changes the topological charge A/| by an integer restricts the 
allowed values of a to the set i?iZ. In terms of their action in the topological charge lattice, 
the vertex operators with a = jR\ for j = 1,2,... (f\ — 1) take us between the different 
sectors of the 4>i field. 

The partition function and correlation functions can readily be obtained since T~Lq is 
quadratic. If we imagine that our system is connected to a reservoir (e.g., a normal metal) 
such that the edge can gain or lose electrons without disturbing the bulk state, then the 
physically meaningful partition function is obtained by summing only over states in a single 
sector. Going between sectors requires the addition of an edge quasiparticle, which must be 
accompanied by the creation of a bulk quasihole, for which there is a finite energy gap. At 
finite length L and inverse temperature j3 the partition function in the ground state sector 
(i.e., the sector containing 1 , 0) ) can be written 

Z ((3) = Tr(e-^o) = ^xT'^ ^) 

+ xg m+n) ' m+n) (gc)xg m - n) ' m - n) (gn), (23) 

where q\ = e~ 2nl3v ^ L , and 

1 00 

X^^^^rE^^ 2 ' with ^( q )^Y[(l- q k ), (24) 
fvl) pe z k=i 

The function x~r^ (q) is the character of a radius R compact chiral boson whose topological 
charge is restricted to the set rZ + s. We note in passing the well known result that Xi (?) 
is equivalent to the character of a chiral Dirac fermion on a circle with antiperiodic boundary 



conditions [i.e., ijj(x + L) = —ijj(x)] Ul7|,|18 



When writing correlation functions we work at zero temperature, in imaginary time 
r = it, and in the limit L — ► 00, a — > 0. Since this involves evaluating ground state 
expectation values of operators that preserve the parity of N c + A^ n , the gluing condition 
does not complicate the calculation. We choose a non-standard engineering dimension for the 
electron operators so that their correlation functions are well behaved in the limit L —>■ 00, 
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%(x) = ^e^W. (25) 



The electron propagator is 

-Qifax) = (T T :*i(T,x): :*J-(0,0):> 

= r~r^ ? r- (26) 

[27r(f c r - zx)]^ - [27r(w n r - ix)]~^~ 

The equal-time Green's function has a spatial dependence Qu(0,x) ~ x~ m which clearly 
signals Luttinger liquid behavior for m ^ 1. The two-point function of the density operator 
©is: 



T>ij(T,x) = (T T pi(T,x)pj(0)) 



1 1 1 

+ oL I\ r^w (27) 



2(m + n) [2tt(v c t - zx)] 2 2(m - n) [2n(v n r - ix)\ 

Note that the off-diagonal term, T>\2 = T>2i, is non-zero even in the absence of tunneling; it 
vanishes only when there are no interlayer correlations in the bulk (n = 0) and no interlayer 
interactions at the edge (g = 0). 

B. Interlayer Tunneling and the Quantum Chiral sine-Gordon Equation 

We now consider adding interlayer electron tunneling to the edge theory. Interlayer 
quasiparticle tunneling at the edge is not considered since the addition of a quasiparticle to 
the edge requires the creation of a bulk quasihole, and hence it is not a low-energy excitation. 
The electron operators defined above have the property that 

e iUi{x) e iUj{x') _ e iUj(x') e iv,i(x) e iw sgn(x' -x)Kij ^ f QJ . % x > ^8) 

If K\2 = n is even, this implies that the electron operators in different layers commute rather 
than anticommute. This can be remedied by the introduction of Klein factors constructed 
from the topological charges in such a way that the tunneling Hamiltonian defined below is 
not modified. The details of this procedure are given in Appendix |A]. 
The Hamiltonian describing interlayer single-electron tunneling is: 



r L / 2 r -, 

Hx = A J L2 dx [ : ^x(x) * 2 f (x) : + h.c.J , (29) 



where A is the tunneling amplitude, which is assumed to be uniform along the edge. Using 
the transformation ( fL5|) and removing the normal ordering with Eq. (|19|) , the full Hamilto- 
nian TC = TCq + 7~Li can be expressed in terms of the rotated Bose fields as 



-L/2 

H = dx 



-L/2 

/3 2 



1 1 2A 

-j-u c ■■ (d x <f) c ) 2 : +~ r v n :{d x (j) n ) 2 : +— — cos(^ n ) 
47T 47r (27ra)^ > 2 



nM+Xsai^], (30) 
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where (3 = \j2(m — n), and A = XoL n has length dimension 2 /2 — 2). Note that because 

of our normalization of the bosonic fields, (3 differs from the inverse period conventionally 
used in the SG theory by a factor of \[2~tx. 

The Hamiltonian separates into two commuting pieces, a free chiral boson Hamiltonian 
(TCf) for the field C , and a chiral sine-Gordon Hamiltonian (H^ SG ) f° r the field 4> n . If we 
had taken v ^ v ' there would be a cross term in the Hamiltonian of the form d x (p c d x (j) Q , 
while if we had taken m ^ ml the two pieces of the Hamiltonian would not commute. For 
the case olv^v' we could perform a further orthogonal transformation on the Bose fields 



which would preserve the commutation relations (16 ) and eliminate the cross term in the 
Hamiltonian, but then the tunneling term would involve both transformed fields M. 

An important property of the tunneling Hamiltonian, Hi, is that it commutes with Ho- 



This follows from Eq. (|C11[) in Appendix |C[ Physically, it is a consequence of the fact that 



the translationally invariant perturbation conserves momentum, which in a chiral system 
with a single velocity is proportional to the energy. It follows that first-order degenerate 
perturbation theory for this system is exact, and therefore the exact eigenvalues of H^sg 
depend linearly on A. 

Since the zero-energy ground state of Hq, |0, 0), is non-degenerate, another consequence 
of the vanishing commutator, [TL ,TLi\ = 0, is that 1 , 0) is also an exact eigenstate of the 
tunneling Hamiltonian. In fact, Hi annihilates the ground state of Hq, as we now show. 
Using the mode expansion ([[7]) we can write 



if3 2 /2 V r»\ /y 2 V L J 

X 



e -i/W^ e *M +) (a:) e Wl ) (x) e 27T0R n N n x/L _|_ _> 



(31) 



where we have introduced 



1 



= —e- k ^ 2 e lk ^b nr , = (0H)t. (32) 



r=l 



If we act on the state |0, 0) with Eq. ( |3TD and imagine expanding the exponentials involving 
N a and (f)t~> we see that only the zeroth-order term survives. We then have 



x 



dx cos 10, 0) = I ^ (~f) e ~ l$XjRn J m dx e ^ X/L 

+ (/3^-/3)l|0,0), (33) 



a, r 



1+ee-e e( 1 ] E ^ Sj — 7 w n b i 

r=l«i=l s r =l ^Vllj=i j=l 



where we have expanded the exponential of 4>^ ■ The a; integral vanishes in every term 
(since (3 2 G Z^ en ), establishing that |0, 0) is an exact zero-energy eigenstate of 7i + 'Hi- Of 
course this does not imply that this state is the ground state of the full Hamiltonian. 

The chiral sine-Gordon Hamiltonian we have arrived at from the bilayer quantum Hall 
system has the properties: /3 2 G Z^ en , and i? n = 1/(3. It is natural to ask if the xSG theory 
can be sensibly defined for more general choices of these parameters. As discussed above, the 
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requirement that cos(/30 n ) changes the winding number by an integer gives the restriction 
(3/ R n = p for some peZ. Next consider the behavior of the perturbation under x — > x + L, 

cos /30 n (x + L) = (e^We^ w + e -tfM*) e -2*tf*«^ _ (34 ) 

Since this operator appears in the Hamiltonian we demand that the periodicity be indepen- 
dent of the state on which it acts, which implies (3R n = p' for some integer p'. Combining 
these two conditions yields (3 2 = pp', and R 2 = p'/p. Without loss of generality we can take 
p and p' to be positive. We find that in addition to the case (3 2 G Z^ en previously consid- 
ered, there is the case (3 2 G Z^ dd . We shall refer to these as the "bosonic" and "fermionic" 
sequences, respectively, since from Eq. (|34] ) the perturbation is periodic for even (3 2 and 
antiperiodic for odd (3 2 . 

To reiterate, for non-integer values of (3 2 there does not exist a choice of the radius R n 
that ensures that the Hamiltonian density changes the topological charge by an integer and 
has a periodicity under x — > x + L that is independent of the state on which it acts. For a 
given value of (3 2 , different choices of the radius do not correspond to distinct theories. If for 
a fixed value of (3 2 we take the minimum value of R n , we recover the condition R n = 1/(3. 
The other allowed radii at this value of (3 2 are integer multiples of this minimum radius and 
by Eq. (|2~3]), they correspond to restricting the allowed values of the topological charge in 
the theory with the minimum R Q . 



III. RG FLOW OF THE CHIRAL SINE-GORDON THEORY 



At the end of the previous section we found that the %SG theory appears to be sensibly 
defined only at discrete values of (3 2 . This suggests that (3 does not flow under a renormal- 
ization group transformation. In this section we examine this question directly. We work in 
the limit L — > oo and consider the Euclidean action of the xSG theory jTD| . 



S 



X SGl 



dx dr 



' 1_ 
Ait L 



id T (j)d x (j)+ (d x <j>)' 



+ 



A 



(2vra) 



cos(/30) 



(35) 



where A = 2A(27ra) 2 ~ /32 / 2 is the dimensionless coupling and in this section we set v n = 1 
for simplicity. Here and in the remainder of the paper we will suppress the subscript on 
the neutral boson, i.e., <fi = <p a . Recall that the action of the standard (i.e., non-chiral) 
sine- Gordon theory is p| : 



Ssg[4>] 



dxdr ( — 



(9 T 0) 2 +(40) 2 i+4 c ° s (^)V 



(36) 



The RG behavior of S$g is the well known Kosterlitz-Thouless flow pi| -|24 

d(3 2 



d\ _ 
d£~\ 



y)A + 0(A 2 ), 



-c A 2 + 0(X d 



(37) 



where the UV momentum cutoff decreases as £ increases, and the value of the constant 
c > is non- universal. In Ref. H, Renn argues that S x sg can be mapped into Ssg by a 
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coordinate transformation, which would imply that the RG flow of the %SG theory is also 



given by Eq. (37). This argument is flawed since at the level of the free action, A = 0, 
the two theories clearly have different numbers of degrees of freedom (the non-chiral boson 
can be written as the sum of two chiral bosons) and hence cannot be related by a change of 
coordinates. The classical (real-time) equation of motion for a free chiral boson which follows 
from the Minkowski version of the quadratic part of Eq. (|35|) is (d x dt + d 2 )<p = 0, which 
has a general solution <p(t, x) = 4>r(x — t) + $(t) where $(t) is a gauge degree of freedom. 
The coordinate transformation Renn proposed amounts to promoting this unphysical degree 
of freedom into a left-moving boson, which when combined with the physical right-moving 
degree of freedom, <f)R, gives a non-chiral boson. 

A standard perturbative, momentum-shell RG analysis for S x sg gives to first order in A 

& 0. (38) 




dt 

We see that tunneling is irrelevant for f3 2 > 4 and relevant for f3 2 < 4. The irrelevant region 
includes all fractional quantum Hall states without interlayer correlations, i.e., n = and 
m — 3, 5, . . . Whether an 0{\) term in the flow of p 2 is present (as it is for the non-chiral 
theory) is a difficult question to answer unambiguously using a Wilsonian RG. Indeed, to 
rigorously obtain the — CqX 2 term in Eq. (|37| ) from a momentum-shell RG one must use a 
smooth cutoff ||2T|. Only because of the Euclidean (rotational) symmetry of the free part of 
Ssg can this term be found simply by using an ad-hoc regularization of the divergent integral 
that would be rendered finite if one employed smooth momentum-space slicing pl^p4|. The 



free chiral boson action does not possess this symmetry, and even continuous regularizations 
lead to different results depending on the details of the particular cutoff procedure. 

Nevertheless, we now argue that there is no flow in (3 for the x^G theory, i.e., dp 2 jdi = 
to all orders in A. Consider calculating the correlation functions of some vertex operator 
perturbatively in A, for example the two-point function: 

^/ :e ia^) ::e -i^).\ : !__J 1 + _^L [d 2 Wl d 2 w 2 1 



X 



- { 

( (gi - wi){z 2 - w 2 ) \ 
\(zi ~w 2 )(z 2 -wij) 



zi - z 2 Y { 8a 4 -/ 32 J (twi - w 2 y 2 



+ (wi <-> w 2 ) - 2 



+ G(A 4 ) , (39) 



where z = r — ix, and /i is an infrared momentum cutoff needed since we are working at 
L — ► oo. We see that for general a the multi-valuedness of the correlation function is altered 
by the perturbative corrections. This should be contrasted with the standard SG theory, 
where (w — z) in Eq. (|39|) is replaced by \w — z\ and correlation functions of all vertex 
operators are single-valued at every order in A for arbitrary f3. Note however that for the 
set of allowed vertex operators defined in Section [II], a — p/(3 where p 6 Z, we see from 
Eq. (^) that the 0(X 2 ) term has the same multi-valuedness as a function of (z — z') as the 
leading-order term. This property holds for iV-point correlation functions (N > 2) and for 
higher-order corrections in A, since the factor of f3 cancels whenever we contract an external 
vertex operator [e 1 ^/^)^)] with an internal one \e l ^ Wi >\. It is conceivable that the theory 
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could be denned at non-integer values of p 1 for this subset of operators, but we still contend 
that there is no flow in (3. 

Suppose we begin with some bare values of the parameters, A and /3 , defined for a theory 
with a UV momentum cutoff A ~ 1/a and imagine running the RG so that the cutoff is 
reduced to A/s where s > 1, and Ao — * A(s), $q — > $(s). We allow for a field renormalization 
by defining (fi'(z') = Z^ 1 (s)(p(z), where z' = z/s. The covariance of correlation functions 
under the RG flow implies that for \z[ — z' 2 \ 3> (A/s) -1 , 



g ia0(*i) . . e -iacf>(z 2 ) A _ /. giaZ^s)^' (z[) . . ^iaZ^s)^ 1 (z' 2 ) 



(40) 

where the l.h.s. is evaluated in the bare theory and the r.h.s. in the renormalized theory. A 
key point is that any flow in f3 is governed by the field renormalization: (3{s) = Z^s)^. If 
we choose a = p/(3q for some integer p, so that the perturbative corrections on the l.h.s. of 
Eq. ( |4Tj| ) do not change the multi-valuedness of the leading-order term, then the requirement 
that the same thing be true on the r.h.s., along with the relations Z^(l) = 1, = $o, 
implies aZ^s) = p/(3(s). Eliminating Z$(s) and a from this equation yields (3(s) = /3 , i.e., 
13 does not flow. This argument for the non-renormalization of (3 can at best be viewed as 
heuristic; its validity is intricately connected to the way the RG transformation is defined. 
We believe that at least the weaker statement that there exists an RG procedure for which 
there is no flow in /3, is true. The exact solutions of the theory at f3 2 = 2 and p 1 = 4 
presented in the next section support this claim. 

In light of the argument that d/3/d£ = 0, the flow in A at the lowest-order marginal 
point, p 1 = 4, will be governed by the higher-order terms in the first expression in Eq. (j38[). 
Because of the complexity of the higher-order bosonic perturbation theory, we shall postpone 
this calculation until we find a fermionic representation of the theory in the next section. 
We will find that at this value of p 1 the tunneling perturbation is exactly marginal, i.e., there 
is a line of fixed points parameterized by A. The RG flows of the standard SG and %SG 
theories are shown in Fig. El. 



X 






t 1 i 




i 1 1 _ 







p 2 







FIG. 3. Left: the perturbative RG flow for the standard SG theory. Both A and p 1 are 
renormalized. Right: the perturbative RG flow for the xSG theory. (3 does not flow and there is a 
line of fixed points at p 2 = 4 (dots). 



IV. EXACT SOLUTIONS OF THE CHIRAL SINE-GORDON THEORY 



In this section we present exact solutions of the %SG theory at the points p 1 = 2 and 



(3 2 = A. 
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A. Relevant Tunneling, p 1 



At the point p 1 = 2, which corresponds to m — n — 1, the cosine term is relevant. There 
is an infinite sequence of K matrices describing maximally chiral edges which satisfy this 
condition: (m, m', n) = (1, 1, 0) + r(2, 2, 2), r = 0, 1, . . . The discussion here also applies to 
the bosonic states (m, m', n) = (2, 2, 1) + r(2, 2, 2). 

At this value of (3 the vSG Hamiltonian is 



njP 2 =2 
n xSG 



L/2 
-L/2 



dx 



1 



A 



in/20 



+ h.c. 



(41) 



where 4>{x) is a radius i? n = \j\p2 chiral boson. A similar problem, with A replaced by a 
Gaussian random variable, appears in Kane, Fisher, and Polchinski's work on the disordered, 
single- layer, v = 2/3 edge 0. A special property of a radius \j\[2 free chiral boson is that 
its character ( f24"D can be expressed as: 



(1/2) 



(42) 



where ^( 2 )i(?) ^ s the character of the spin-j irreducible representation of the su(2) 1 Kac- 
Moody (KM) algebra EH]. Indeed, the trio of dimension one operators d x (f>, e 1 ^^, e~ lv ^^ are 



related to the su(2) 1 generators, J a (x), a = x,y } z } by |2 

1 



J z (x) 



-.d T d(x) 



2ttV2 

In terms of the Fourier components 

r L/2 



j±( x ) = r±ijy 



2na 



D ±iV2<t>{x) 



J"= / dxJ a (x)e 

J-L/2 



-2-Kinx/L 



71 G 



the algebra is 



1° ni u n 



n 



5 ab 5, 



n+mfi 



+ ie abc r 



n+rn' 



With the help of Eq. (pE3|) the Hamiltonian (II) can be expressed in terms of currents 

rL/2 



n 



/3 2 =2 
X SG 



L/2 



dx {2ttv q :[J z {x)f: +2\J x (x)} 



We now use an identity valid for any state I7) in the Hilbert space 

r L / 2 9 r L / 2 1 o 

/ dx :[J z (x)] 2 : 7 = / dx- :[J(x)] 2 : 7 , 

J-L/2 J-L/2 3 

to write Eq. (f46"|) as: 



n 



2 =2 
X SG 



L/2 dx (2^5 . [J(;r) ]2 . +2A J*(x; 

-L/2 13 



(43) 



(44) 



(45) 



(46) 



(47) 



(48) 
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One way to prove Eq. ([47]) is to note that the operators J z obey a u(l) subalgebra of the 
su(2)i algebra generated by the J% operators. We can construct two families of operators 



T = — ■ T a T a 

~ 3 ' ~ 



r (z) _ . jz j 



m"n—m " ' 



(49) 



each of which obey c = 1 Virasoro algebras. Since the J z generate a subalgebra we have 
[L n , L)$\ = [L^\ L$], and therefore by the standard Goddard-Kent-Olive (GKO) coset 
construction ^6] the operators L$ = L n — L)fi> obey a c = Virasoro algebra. It is well 
known that the only unitary representation of the c = Virasoro algebra is trivial, i.e., 

It) = f° r an n P^C - We thus conclude that LJ7) = I7) for any state in the Hilbert 
space, and for n = we obtain Eq. (|47D . 

Given any matrix R G 50(3), we see from Eq. ( f4"5p that the currents J a (x) = R ab J b (x) 
also obey an sm(2)i KM algebra since 5 ab and e abc are S0(3)-invariant tensors. If we make 
the particular choice R ab = S ay 6 by + S az S bx — S ax S bz , which corresponds to a rotation by 7r/2 
about the y-axis, and express the rotated su(2)x currents in terms of a new radius R = l/y/2 
chiral boson, 9(x) (i.e., J z (x) = [d x 8(x)]/(27T\/2), etc.) then from Eqns. ( f47l) and (|48|) we 
have 



L/2 

-L/2 
L/2 

-L/2 



eke < 2nVr, 



J z {x) : +2AJ z (x) 



dx 



^v n :(d x 9f:+^d x 9 

47T Z7T 



(50) 



The mapping between Eqns. (p|) and (0), when expressed in terms of the two Bose 
fields corresponds to the non-linear transformation 



d x 6(x) = — cos (y/2(p(x)), d x <j)(x) 
a 



a/2 

-— cos(v / 2^(a;)). 
a 



(51) 



While we could have taken the first expression in this equation as our definition of 8(x) and 
proceeded by showing this implies : (d x <p) 2 : = : (d x 8) 2 :, the detour through the KM algebras 
is more deductive. 

We have succeeded in writing the Hamiltonian as a quadratic form in terms of a single 
chiral boson. The tunneling term is linear in the topological charge of the field 9(x), and 
thus it is clear that the A = ground state, which is annihilated by the perturbation, is 
not the ground state for A / 0. However, since the A = ground state remains an exact 
eigenstate of Eq. (|50D for all A [see Eq. (|3~3D1 we conclude that zero-temperature perturbation 
theory in A gives incorrect results. In terms of first-order degenerate perturbation theory for 

the low- lying eigenvalues of 'H^sg '■> this is a consequence of a level crossing for A = 0(1/ L). 

Several authors have claimed that a gap is generated when the perturbation is relevant 
|i~l]2"8|j. It is clear from Eq. ([50]) that this is not the case. We still have a (gapless) chiral 
boson, whose topological charge sectors are reorganized in energy. The energy spectrum of 
the winding number mode for two different values of A is shown in Fig. [|. For the points 
A = (jrv n /L)M, where M G Z, this reorganization can be reduced to an overall shift in the 
energy spectrum by defining Ng = Nq + M. The subtlety is that while we have a gapless 
mode for A = (4>(x)) and a gapless mode for A 7^ (9(x)), the modes are non-trivially 
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connected even for arbitrarily small A. We see from the relation, Eq. fl5T|) , that exciting 
modes of the original boson involves the operator d x <p ~ cos(v / 2^), which changes Ng and 
hence involves an energy of order A even for small momentum excitations. Therefore the 
original field appears massive, whereas the number of gapless modes is unchanged. 



\ \ 


E N 

/ 

/ 

// 






N 


I i 







FIG. 4. Energy spectrum for the topological charge (N) of the field 9{x) at two different values 
of A. The open circles indicate the ground state. For small A (solid line) the ground state is 
unchanged, while at a larger value of A (dashed line) the ground state has a non-zero winding 
number. 



From Eq. (|50| ) we can immediately write down the partition function for the x^G theory 
in isolation 

Z" x 2 sc 2 (P) ~E e-^ q f/\ (52) 

To construct the partition function of the full edge theory we need to consider gluing con- 
ditions. Recall from Section [II] that we need to separate states according to their A^ n -parity, 
as in the partition function at A = (p3|). For A ^ 0, eigenstates of the Hamiltonian are 



no longer eigenstates of N n . However, we see from Eqns. ((43|) and (f44"D , that N n = 2Jq and 
similarly Ng = 2Jq, where Ng is the winding number operator of the field 8(x). Therefore, 
states in the spin-0 (1/2) representation of the J% algebra have N n even (odd). Under the 



rotation from J a to J a the two representations do not mix, and hence eigenstates of Eq. ( p0|) 
which have Ng even (odd) are linear combinations of iV n eigenstates, all of which have even 
(odd) N n eigenvalues. Thus we can write the partition function including the charged mode 
as: 

Z $2=2 (P) = xT +n)fl \^s = c ' 6Ven (/5) + xT +n) ' m+n \^s = c '° dd (/3), (53) 
where Z^ 5 Q 2 ' even is given by Eq. fl52"D with the restriction N e Z even , and similarly for 

7 /3 2 =2,odd 
Z XSG 

We now turn to the computation of correlation functions. We define the tunneling 
operator as the Hamiltonian density of the tunneling perturbation 

G(x) = srr cos \3d>(xj\ , (54) 
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and denote its two-point correlation function by T(r,x) = (T T 6(r, x)0(O)). Using transfor- 
mations (|l^) and (pi]) we can express the charge density (|TTD , and tunneling operator ([54|), 
in terms of the fields C and 

I i 

piA x ) = i ^d x (j) c (x) ± - — cos V29(x), (55) 

/87r 2 (2m - 1) 2 ™ 



0(x) = -^^(x). (56) 

7TV 2 

Since the Hamiltonian 7Yf[0c] + ^xSG 2 ^] i s quadratic in these fields we find 

V-(r x) = 1 1 + cos ( 2Xx / v n) (57) 

13 K ' ; 2(2m - 1) [2tt(v c t - ix)f 2 [27r(w n r - ix)f ' V ' 

T{t,x) = o + ■ (58) 

The part of the density-density correlation function arising from the neutral mode is spa- 
tially modulated when tunneling is present, while the two-point correlation function of the 
tunneling operator acquires a constant piece. Both of these are a result of the non-zero 
topological charge of the 9 field in the ground state. 

The calculation of the electron propagator is slightly more involved. From Eqns. ( |T5] ) 
and (|25|) we find that for m — n = 1 the electron operators can be written 

= ■jk^^"^*' '• (59) 

and therefore they cannot be readily expressed in terms of <p c and 9. However, from this 
expression it is clear that the single-electron Green's function factorizes: 



n F [4>c] 

xje^H"^) ^ . (60) 



H xSG 



The correlation function involving the <fi field is the same for all values of m and n that 
satisfy m — n — 1. The correlation function involving </> c can be computed for all m and n 
since 7"^f[0c] is quadratic. 

For the special case m = l,n = the l.h.s. of Eq. ( |60|) can be calculated exactly as we 
now demonstrate [p9fl . For the uncorrelated integer 110 state there exists a chiral fermion 
description of the edge theory including tunneling, TC + Ti.%: 

H {110) = [ dx : [-w^Ja^ + 2tt#I'0i'02^2 - A(V»JVi + h - c -)l : • (61) 

J-L/2 L J 

The details in going from the bosonic form of the Hamiltonian ([3(]) to the above fermionic 
form are discussed in Section |V|. If we define ip±(x) = (ipi(x) ±ip2(x)) /y/2, Eq. fl5"T|) becomes 
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L /2 

H {m = [ L dx: [-iv (jp\d x ip + + ipld x ip-) + 2ng^ + ip+^_^ - \{ip\.'ip+ - iplip. 



(62) 



This Hamiltonian can be brought to a quadratic form by bosonizing according to ip±{x) = 
e «/>±0) l^/2ixa. The resulting quadratic form is diagonalized by defining 9 + = ((f> + +0_)/-\/2, 
and 9- = — <p_)/y/2 + \[2\xjv a . The exact single-electron Green's function is then 
readily found: 

(no), , _ 5jj cos(\x/v n ) - (1 - Jjj)isin(Ax/u n ) , . 

— y i7 - (r, xj — / • K 106 ) 

2ti{v c t — ix)2ir{v n T — ix) 

We can thus use Eqns. (|60|) and ( |63f) to find the single-electron Green's function for all states 
in the m — n = 1 sequence: 

SijCosjXx/vn) - (1 - ^j)zsin(A3;/^n) 

yij{T,X) — l ; — ) lo 4 J 

[27r(t> c r — zx)J 2 J2^{v n T — ix) 

the only modification being the exponent of the piece corresponding to the charged boson. 



B. Marginal Tunneling, p 1 = 4 

The point (3 2 = 4, at which the cosine term is marginal, corresponds to the case m— n = 2. 
There is again an infinite sequence of K matrices describing maximally chiral edges which 
satisfy this condition: (m,m',n) = (3,3,1) + r(2,2,2), r = 0,1,... The results here also 
apply to the bosonic states (m, m', n) = (2, 2, 0) + r(2, 2, 2). 

At this value of (3 the %SG Hamiltonian is 



n 



X SG 



L/2 
-L/2 



dx 



47T 



A 



(2na) 



(e 2 ^ + h.c. 



(65) 



where <f)(x) is a radius R n = 1/2 chiral boson. If we restrict ourselves to the sector of the 
edge theory that contains the ground state, then we find that N n G Z even , since adding an 
electron to layer 1(2) changes N n by =p2. From the expression for the compact chiral boson 
character (|2"4[), it is clear that a radius 1/2 boson restricted to even topological charges is 
equivalent to a radius 1 boson with no restriction [i.e., Xi/'^il) = Xi* (<?)]• Hence, for this 
sector of the edge theory we can take <p(x) to have radius 1 and thus define a chiral Dirac 
fermion ip(x) = e l ^ x ^ / y/2na. This result, that the edge theory of double-layer states with 
m — n = 2 contains a chiral Dirac fermion was previously established in the wavefunction 



approach by Milovanovic and Read [14[ . 

To find a fermionic representation of the tunneling term, we first consider bringing to- 
gether two chiral Dirac annihilation operators 



ip(x)ijj(x + e) = eip{x)d x ^{x) + 0(e 2 



(66) 
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where the zeroth-order term vanishes by Fermi statistics. We next consider bringing together 
the bosonized representation of the Fermi fields 



2ixa 



L 2 
2-Kie 



: e 



2i(j>{x) . 



+ 



L 2 



:e 2 ^: -a 



o x (p{x) + — 



:\ + 0(e\a 2 ). (67) 



Since the terms on the r.h.s. are normal ordered, we can safely take the limit a — > 0, and 
compare the order e terms in Eqns. (|66|) and (|67|) to obtain the identity 



if)(x)d x i(>(x) 



2na 2 



(68) 



We verified the validity of this formula by checking that the iV-point correlation functions 
of the operator on the l.h.s. (and its Hermitian conjugate) calculated in the free Fermi 
theory are identical to the Appoint correlation functions of the operator on the r.h.s. (and 
its Hermitian conjugate) calculated in the free Bose theory. The proof of this equivalence 
is given in Appendix [FJ. In Appendix |C| we demonstrate that the correspondence holds 
for the commutation relations as well. Readers untroubled by operator manipulations in 
bosonization can ignore them. 

Using Eq. (|68|) we can then fermionize the x^G Hamiltonian fl65D to 



n xSG 



L/2 
-L/2 



dx 



-iv n ip^d x ip - i^-^d^ + ipd x ip) 



(69) 



The effect of the tunneling term is made transparent by writing the Dirac fermion in terms 
of its Majorana components: ip{x) = [xi(x) + 1x2(2)]/ V2, where x\ = Xu 



H 



1 r L / 2 



X SG 



2 J-L/2 



dx 



(70) 



An alternative derivation of this form of the Hamiltonian is given in Appendix 0. We see 
that the effect of tunneling is to split the neutral mode into two Majorana fields with different 
velocities. This is a c = 1 = (1/2) + (1/2) analog of the usual spin-charge separation in one 
dimension [c = 2 = (1) + (1)]. Since the Hamiltonian (|69"D is quadratic, modes of the fermion 
field with momentum \k\ > A, for some A, do not couple to modes with momentum \k\ < A, 
and thus the tree- level RG flow is exact. Therefore if we define our RG transformation so 
that the A = Hamiltonian is a fixed point, we immediately see that the coupling A is 
exactly marginal, and we have a line of fixed points as asserted at the end of Section |ITT[ 

Armed with a quadratic representation of the Hamiltonian we can readily compute the 
partition function and correlation functions exactly. From the expression for the partition 
function (p3|) at A = we see that the characters of the neutral mode that appear at the 
(3 2 = 4 point are Xi$ an d X1/2 ■ These can be written in terms of the characters of the 
highest weight h representation of the Virasoro algebra with central charge c, denoted V£: 



(4 '%) 



Xi/2 



V 1/2 (<?)] +[v 1 1 / / 2 2 (g) 



(4 < 2) (?) 



Xi/2 



2VQ 1 



1/2 ( ? )V^ 2 (g). 



(71) 
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Since a single chiral Majorana fermion with antiperiodic boundary conditions has a character 
Vq^ 2 + Vl/2i we see fr° m Eqns. fl23|) , (|70f) , and (|7lD , that the partition function of the %SG 
theory in isolation is 



V 1/2 (g + ) + V$(g + )] [v 1/2 (g_) + V$(g_ 



(72) 



where g± = e 2 ^P V ±/ L ^ and = v n ± A/7T, while the partition function of the full edge 
theory with tunneling, including the charged mode with the appropriate gluing condition, is 



^ 2 = 4 (/?)= X g m+n) ' 0) ' 



(9c) [Vo /2 (q + )Vo ,z (q-) + V$(g+)V$(g_) 



,1/2 



1/2, 



;V2, 



(2(m+n),m+n) , 



(?c) [V 1/2 (g + )V^(g„) + V^(g + )V i/z (g-) 



,1/2 



1/2, 



(73) 



Recall that there are three irreducible Verma modules in the c = 1/2 minimal model 
with highest weights h = 0, 1/2, and 1/16, which in the Ising model terminology correspond 
to the identity (I), energy (e) and spin (a) fields. The characters of the first two of these 
representations appear in the above partition function, the remaining one occurs if one 
considers the other sector of the neutral mode which involves the characters X1/2 an< ^ 

X1/2 j an d which corresponds to the addition of a quasiparticle. These combine to give 

(4,1) (4,3) _ 2 / v l/2\ 2 

Al/2 + Al/2 — Z \ V 1/Wj ■ 

Turning now to correlation functions, we can use the transformation flT5|), and the 
fermionization of the neutral boson to write the electron (25), density (|TID , and tunnel- 
ing (K4T) operators in terms of the fields C and ijj, 



Pi 

e(x) 



1 



(27ra)(™- 1 )/ 2 
1 



d x 4> c + 



-1 



i 

~2^ 



6 a ^ + 5 i2 ip 



: (V ^ + ^d x if) : , 
where -z = 1,2 is the layer index. The single-electron Green's function is 



- Qij(r,x) 



2 [2tt(v c t — ix)] 



m—l 



+ 



2n{y + T — ix) 2tt(v-T — ix) 



(74) 
(75) 
(76) 

(77) 



and we see that in the presence of tunneling an electron has broken up into three pieces, 
propagating with three different velocities. The two-point functions of the density and 
tunneling operators are: 



T>ij(r,x) 



T(t'x) 



1 



4(m - 1) [2tt(v c t - ix)} 
1 

+ 



2 + 
1 



■!)*+> 



1 



4 2ii(v + t — ix) 2tt(v_t — ix) ' 



[27t(v + t — ix)] [2tt(v-T — ix)] 



4 • 



(78) 
(79) 
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V. FERMIONIZATIONS OF THE CHIRAL SINE-GORDON THEORY 



In this section we describe a general method for fermionizing the xSG theory for j3 2 G Z + 
by the introduction of auxiliary degrees of freedom. The procedure is illustrated in detail 
for the cases (3 2 = 2 and p 1 = 4. We also show that for these two points the fermionized 
theory with auxiliary degrees of freedom can also be solved exactly, and the projection onto 
the xSG Hilbert space can be performed explicitly reproducing the results of the exact 
solutions given above. In terms of the solution of the double-layer problem, the method 
described here is redundant. However, we believe it may be useful in the study of quantum 
Hall systems with more than two layers. This section can be skipped since the subsequent 
sections are independent of the material presented here. 



A. General Strategy 



In the case of the standard SG theory it is possible to fermionize the theory at all 
values of p 2 < 4 to the massive Thirring model [pOfl . For the %SG theory we have found 
fermionizations at the discrete points (3 2 G Z + . Recall that all the points in the bosonic 
sequence (/3 2 G Z+^J correspond to some double-layer quantum Hall system, while none of 
the points in the fermionic sequence ((3 2 G Z^~ dd ) do. The scaling dimension of the tunneling 
perturbation is an integer in the bosonic sequence and a half-integer in the fermionic one. 
We defer discussion of the fermionic sequence until Section |VT[ 

The algebraic identity at the heart of fermionization (or bosonization) is the Cauchy 
formula for the determinant of an N x N matrix M i? - = l/(zi — z' rj ) which is a function of 
2N complex variables z^z'^i = 1 



-l)N(N-l)/2 detM (N) 



lll<i<j<N{ z i — z j){ z 'i — z j) 
Y[l<i,j<N\ Z i ~ z j) 



(80) 



The determinant on the l.h.s. is proportional to the 2iV-point correlation function of a single 
chiral Dirac fermion: 



N 



W=l 



N \ f_lW(^-l)/2 



(81) 



where Zi = (vr^ — iXi), and similarly for the primed coordinates, v is the velocity of the 
fermion, and time ordering will be implicit from here on. We desire a fermionic representation 
of the tunneling operator, Eq. (0). The 2iV-point correlation function of the tunneling 
operator is a linear combination of terms of the form 




-iP<t>(T!, x r) 



ni<i<j<N( z i ~ z j)( z 'i ~ z j) 



Hl<i,j<N\ z i z \ 



detAfW 



12) 



For the bosonic sequence it is clear from Eq. (|81 
be reproduced by an identification of the form 



that the (3 2 power of the determinant can 
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JP<t>(x) 



/? 2 /2 



n m*) 



B 2 

n 



ip\(x), 



(83) 



where the ipi(x) are p 2 independent chiral Dirac fermions with identical velocities. This 
construction closely resembles Jain's bulk parton construction for hierarchical fractional 
quantum Hall states pO and its edge analog |12|, although the motivation here, a desire for 



a fermionic representation of the tunneling operator, is very different. 

In the xSG Hamiltonian there is only a single chiral boson (with central charge c = 1), 
so we first have to add (f3 2 — 1) auxiliary c = 1 degrees of freedom (e.g., chiral bosons) in 
order to carry out this procedure. Of the points in the bosonic sequence, only the first two 
(f3 2 = 2,4) occur in the region where the cosine term is not irrelevant [by Eq. flSsDI, and 
these are the ones we now describe in detail. 



B. Relevant Tunneling, p 1 



As noted above, an additional degree of freedom is necessary for a straightforward 
fermionization of the x^G Hamiltonian, so we add to the Hamiltonian an auxiliary free 
compact chiral boson ((f)) with radius R, topological charge N, and a velocity equal to the 
velocity of 0. The full Hamiltonian is then: 



n 



f3 2 =2 
XSG 



n 



/3 2 =2r 
X SG I 

L/2 



dx 



L/2 



+ 



A 



-U n :(d x (j)f: +^-v Q :(d x 4>) 2 : , . 
Air 47r 2ira 



iV2<j> 



+ h.c.) 



(84) 



We then perform a canonical transformation which mixes the field appearing in the tunneling 
term with the auxiliary boson 



J_/l 1 

V2U -i 



0i 



in terms of which we find 

-L/2 



n X SG 



-L/2 



dx \ — v n : d x 6id x 6i : H 

I 4vr * * * * 2 7ra 



3 i(0i-0 a ) 



+ h.C. 



The topological charge dependent part of the transformation in Eq. (|85D is 
RN 



i=(A« + AM), i=iV n = ±=(N«Rl - N*R% 



(85) 



(86) 



(87) 



where Nf , R\ are the topological charge and radius of the field 9i(x), respectively. We choose 
R = l/v2 so that the requirement Nf G Z is consistent with setting R\ = 1. These choices 
also impose the gluing condition N + N n 6 Z even . With these values of the radii for the 6* 
fields, we can define two chiral fermion operators 



ipi(x) 



V2 



na 



(88) 
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These two flavors of fermions commute, but once again this can be fixed without modifying 
what follows, and the details are given in Appendix |A]. Using this definition in Eq. (j8q ) we 
have a fermionized Hamiltonian 



KySG = dx '■ \-iv n if>ld x if)i - A(V4^i + h.c.)l : • (89) 

J-L/2 



As noted in Section [D], the fermion fields introduced in Eq. fl88|) obey antiperiodic bound- 
ary conditions and therefore the allowed values of the momentum k for the Fourier mode 
operators Cki of ipi{ x ) are k G (27r/L)(Z + 1/2). The fermionic normal ordering introduced 
in Eq. (|89|) is defined with respect to the (non-degenerate) ground state of this Hamilto- 
nian at A = which is annihilated by for k > and c ki for k < 0. By introducing 
ip± = (ipi ± ^2) /a/2 the Hamiltonian is diagonalized: 

L /2 

H^sc 2 = J_ l ^ dx : [-iv n (^+8^+ + ipld^.) - A(^lV+ + ^-)] ■ ■ (90) 

Note that in the case of the 110 state the charged boson (</> c ) has radius R c = 1/V2 
and hence the fermionization can be accomplished without adding an auxiliary degree of 
freedom by replacing <fr by <j) c in Eq. ([85]). The result of this procedure is nothing more than 
the chiral fermion description of the original Hamiltonian of two uncorrelated, interacting, 
v = 1 edges, given above in Eq. fl6T|). For the 110 state the advantage of using an auxiliary 
boson rather than the charged mode in the fermionization is that in the former case the 
theory is quadratic (|90|), whereas in the latter case it is not ([61]) . 

Having solved the theory with the auxiliary boson, we now show at the level of the 
partition function that projecting out the unphysical degree of freedom from the spectrum 



of Eq. (pOl) reproduces the result from the direct solution in Section [TV|. From Eq. ( p0|) the 
partition functions of the x^G theory with the auxiliary boson is 

AT + ,7V_eZ 

To write the partition function including the charged mode we must again consider gluing 
conditions. We must separate states according to their iV n -parity, which because of the 
condition N + N n e Z CV en, is equivalent in their A-parity. Since N is just the total fermion 
number measured with respect to the A = ground state, we see the partition function of 
the full edge theory is: 

Z" 2=2 ((3) = xg m+n) ' 0) (gc)ifc 2 ' even + xT +n) ' m+n \^st° d \ (92) 

where Z^ s ^' even is given by Eq. ( |9"TD with the restriction iV + + iV_ E Z eve n and similarly for 

Z^ s ~^' odd . To project out the radius R = l/y/2 auxiliary boson we added (<fi) in Eq. (|84|) we 
must evaluate 

7/3 2 =2. even yf3 2 =2,odd 

XSG { ) " v (2 < 0) (a ) To ) ( 3 
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for the %SG theory in isolation and 



j_ (2(m+n),m+n) 



7/3 2 =2, odd 
J XSG 



(94) 



for the full edge theory. It is useful at this point to recall the definitions of the Jacobi theta 
functions: 



(95) 



If we define v = i/3\/2n and r = i/3v n /L, we see from the definition of the chiral boson 
character (El) and Eqns. (I9H) that 



e 2 


V 


r) 


V 2iriur inrr 2 








rez+i 


e 3 


y 


r) 


— 2-^i > 










e 4 


y 


r) 













1 „ > (2,1) . - I 



e 2 (o|2f). 



(96) 



Next note that we can rewrite the sums appearing in Z^ SG 2 ' even and Z^ S Q ,odd using 



E 1 

^ 2 



E 

A r + +Af_eZeve, 

E = E ■ 

JV++iV_eZ odd N+.N-eZ 



1 + (_i)^v + +iv_ 
1 _ (_i)^ + +^v- 



(97) 



and hence 



7/3 2 =2,even 
J XSG 

A/3 2 =2,odd 



— {^i^ + ie^)] 2 }, 



(98) 



Finally, using the standard doubling identities [ CT 



6 3 (2i/|2r) 



6 2 (2^|2r) 



[e 3 (Hr)] 2 + [e 4 (Hr)f 

26 3 (0|2r) 
[© 3 (Hr)] 2 - [6 4 (Hr)f 



26 2 (0|2r) 

and Eqns. © and @§) in © and (|§ we obtain 
1 



(99) 



<p(<h) 



[9 3 (2i/|2r) + e 2 (2i/|2f)] 



Z" 2=2 ((3) = x£ m+n) '%c)^e 3 (2P|2f) + x ^' m+n \ qc )-^Q 2 (2u\2f), (100) 
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which by Eq. (|95|) are in precise agreement with the results of the exact solution in the 
previous section, Eqns. ([52]) and (|53[) . We have also checked that correlation functions 
calculated with the fermionic Hamiltonian (|90|) are identical to those found in Eqns. ( |57|) 
and (|5|). 



C. Marginal Tunneling, p 2 = 4 

We next consider the generic multi-flavor fermionization of the xSG Hamiltonian at 



(3 2 = 4. We add to 7~ttsG the Hamiltonian for three free chiral bosons 
and topological charges Nf 



with radii FL 



n 



XSG 



n 



X SG 

I 



+ [01,02,0a] 
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fL/2 f 1 r 

/ dx I — : u n (9. 

J-L/2 47T 



+ ( v * + 7T- J [(^^i) 2 + 
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7T 



2tt / 
2A 



(2vra) 



■ cos[20(x)] 



(101) 



The velocities of the auxiliary fields are completely at our discretion, and the motivation for 
the specific choice made here will become clear below. We perform an orthogonal (canonical) 
transformation which mixes the physical mode ((f)) with the auxiliary bosons: 



(102) 
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In terms of these new bosons our Hamiltonian is: 

c£/2 



n xSG 



L/2 
X 



dx 



2vr/ (2?ra) 2 



^(^1-^2+613-04) 



+ h.c. 



The topological charge dependent part of the transformation in Eq. (|102|) is 



(103) 



±N n = \{R{Nl - B°N° + R°N! - i2jJV}), = i(i2?JVf + 

44 = i(i2?iVf - - iigiV| + R{Nf), R 2 N 2 = -1=(B£N$ + R\N\ 



V2' 



"All 



(104) 



where N?,Rf are the topological charge and radius of the field Qi(x), respectively. We 
again want to use our freedom to choose the quantities Ri to ensure R\ = 1. Clearly 
Ri = R 2 = 1/V2, R3 = 1/2 accomplishes this. Since the 0j fields have unit radius we can 
introduce a quartet of chiral fermion operators 



ipi{x) 



V2 



na 



(105) 
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Once again these fermion operators are defined here without the Klein factors necessary to 
ensure the proper anticommutation relations. The fact that this can be remedied without 
modifying the form of the Hamiltonian is demonstrated in Appendix HI 
In terms of these fermions Eq. (|103|) reads: 



n xSG 



dx : \-iv'J)\d x i)i + ^- Ui^^I + V^lV^i 

-L/2 I 1 1 



(106) 



where v' n = v n + A/27T, and A' = 2A. The fermion normal ordering in Eq. ( |106| ) is again 
defined with respect to the A = ground state which satisfies Cfcj|f2) = c^_ ki \Q) = for 
k > and all i, where k e (27r/L)(Z + 1/2) and Cki are the Fourier components of i/ji(x). 
At this value of p 1 the A = ground state is the exact ground state of the Hamiltonian for 
A 7^ (in contrast to the p 1 = 2 case where the A = ground state is an exact eigenstate 
for A ^ 0, but not the ground state). Assembling our four spinless fermions into a pair of 
(pseudo-)spin-l/2 fermions according to 



"02 



V>4 
^3 



(107) 



the Hamiltonian ( |106| ) can be written as a quadratic form 



H 



X SG 



L/2 



L/2 



dx 



in terms of the currents 



7 ^(J 2 1 +J 2 2 ) + 2 -^(J 2 1 + 3l) + X'JvJ2 



(108) 



Ji(x) =:$L(a;) $i 



(x):, = - :$LW^^W 



(109) 



where a a ,a = x,y,z, are the Pauli matrices. The reason for our choice of velocities in 
Eq. (|101|) is now apparent: along with the transformation in Eq. ( 102 ) it produces a com- 
bination of four-Fermi couplings in Eq. (|106| ) that can be written as a scalar product of 
currents. This is not a unique construction, we can obtain the same final form (|108|) with a 
different choice of velocities by modifying the transformation ( |102| ). 
Defining the Fourier components of the currents 



Ji. 



the algebra is 



L/2 
-L/2 



dxJ,(x)e- 2winx/L , 



L/2 
-L/2 



dxiAx)e~ 2mnx/L , ne 



(110) 



" i,ni j,m 



J T a 



n 



-5ij5 ab 5 n+m fl + i5ije abc J^ n+m . 



Ill) 



We therefore have two u(l) ( Jj) and two su(2)\ (Jj) Kac-Moody algebras. From Eqns. (|109|) 
and ( |110| ) we find that the ground state satisfies 
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Ji, n \n) = J£ n \n) = 0, for 2 = 1,2; 



x, y, z, 



n > 0, 



(112) 



and therefore normal ordering for the current operator modes is defined by moving any 
current with momentum index n < to the right past any current with index m > with 
which it has a non-trivial commutator. 

To solve Eq. ( |108| ), the spin part of which is the Sugawara form of two coupled su(2)i 
Wess-Zumino-Witten models [E5], we first note that we can define: 



= - • J J- 



j — _ . ja ja 

J -'n — g ' " i,m " i,n—m • 



;ii3) 



These four sets of operators obey four independent c = 1 Virasoro algebras. The operators 
J + = J x + J 2 obey an su{2) 2 diagonal subalgebra of su{2) 1 © su(2) 1 , and the bilinear 
operators formed from them 



T (+) — ± . ja ja 



(114) 



obey a c = 3/2 Virasoro algebra. From the GKO coset construction, the operators L^/ 2 ^ = 
LW + - LM form a c = 1/2 Virasoro algebra which is independent of the one formed 
from the J +jn currents, i.e., [L^ 2 \ L^} = 0. Using these observations and the trivial 
identity 



: J Ax) ■ J 2 (x) 



l :[J + (^)] 2 : 3 i :[3 1 (x)} 2 : + :[3 2 (x)] 2 : :[J + (x)] 2 : 
2 4 2 1 3 4 



(115) 



the Hamiltonian becomes 

7TV' 



n xSG 



L 



(2) 



+ 



2ttv' X' 

H 

L 2L 



2K 



3V 
2L 



(1/2) 



(116) 



By transforming from Lffl to -^^ 1//2 ^ we have succeeded in diagonalizing the Hamil- 
tonian. From Eqns. ( |112| ), Q113 ), ( 114 ), and the definition of L^ 2 ^ one finds that the ground 
state is annihilated by C$, L^~\ L^ 2 ^ for i = 1,2 and n > 0. Excited states are obtained 
by acting on \Q) with combinations of Virasoro generators with n < 0. 

We now project out the auxiliary degrees of freedom we added in Eq. ( |101| ). The Hilbert 

space of the %SG Hamiltonian is the subspace of the Hilbert space of ii^SG f° r which the 
auxiliary chiral bosons are in their ground state. We define a state I7) as physical (i.e., 
belonging to the Hilbert space of the Hamiltonian) if it satisfies 



L/2 
-L/2 



dx :{d x <Pi) 2 : | 7 ) = 0, for z = l,2,3. 



117) 



Using Eqns. (|102[) , ( |105| ), ( |107| ), and ( |109| ), these conditions can be written in terms of 
currents as 



L/2 

dx : [Ji(x)\ : 7) 

-L/2 

L/2 2 

dx : JUx) : 7) 

-L/2 1 J 



for i 
0. 



1,2, 



(118) 
(119) 
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We note that these quadratic conditions can be recast as linear constraints. For example, 
using Eq. ( [1101 ) in Eq. ( 118 ) and acting on the left with (7I we find 



i,0 



(120) 



n=l 



where we have used (Ji, n ) — Ji-n- Since the norm is positive semi-definite we see that 
Eq. (|118|) is equivalent to the set of linear conditions Ji >n \l) — for n > 0. The condition 
in Eq. ( |118| ) states that we should project out the u(l) algebras, i.e., acting on any physical 
state with £_„ for i = 1,2, and n > 0, necessarily produces an unphysical state. From the 
Hamiltonian ( |116| ), and Eq. ( |118j ) we have 



«fc 4 |7> = (| 



A' 



<+£l4 +) + 



4n 



3A'\ (1/2) 



|7>, 
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if I7) is physical. The second projection condition, Eq. ( 119| ), is slightly more subtle. We 
are aided by the observation that the operators J* obey a u(l) subalgebra of the su(2) 2 
algebra generated by J+ n - Therefore 



r( 2 ) = _ ■ j z j z 

n 2 ' +,m +,n— m " > 



(122) 



obey a c = 1 Virasoro algebra and we can once again use the GKO construction to define 
operators L'^ 1 ^ = — Lffl which are independent of the L$ (i.e., [L'^ l / 2 \ ]J$\ = 0) and 

•W 



which obey a c = 1/2 Virasoro algebra. Using L Q 
and (11221) we arrive at 



(+) 



L' (1/2) +L ( Z> and Eqns. O, (pll) 



2^ 




(1/2) 



+ vL 



3A'\ (i /2) 



It), 



(123) 



if 1 7) is a physical state. One can readily show [L'^ 1 ' 2 ' , L^ 2 '] = and therefore from 
Eq. (|123|) we see that the xSG Hamiltonian can be written as the sum of two independent 
c = 1/2 pieces with velocities (v' n +X' /4n) = v n +\/n = v + and (v' n — 3A'/47r) = v n — \/n = u_, 
in precise agreement with the expression ( fTOD obtained by the direct solution in Section [IV], 
since a Majorana fermion has a central charge of 1/2. 



VI. FERMIONIC SEQUENCE 

In Section [V] we discussed the generic fermionization of the %SG theory for points in the 
bosonic sequence ((3 2 G 2^ en ), which involves the introduction of (fi 2 — 1) auxiliary chiral 
bosons. This procedure can also be applied to the fermionic sequence (/3 2 G T^ AA ), but the 
subtlety here is that the perturbation is antiperiodic, from Eq. (|34]) we have 

cos (p(j){x + L))=- cos (p<p(x)) , (124) 

which is consistent with the fact that its scaling dimension is half-integer and hence in the 
multi-flavor fermionization it is represented as a product of an odd number of fermions. One 
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can modify the definition of the theory at these points by adding a factor of e~ lnx / L to the 
tunneling amplitude A to make the perturbation periodic. For the case (3 2 = 1 the chiral 
boson has radius R n = 1 and the theory can be fermionized without adding auxiliary degrees 
of freedom. With the modification to the tunneling amplitude one finds 
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dx 
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: +Av / L(ci/ 2 + c\ 



/2J 



(125) 



reZ+l/2 



where we have used e l ^ x > / y/ 2ira = ip(x) = L^ 1 l' ls £ JT e lmrx l L c r . A similar Hamiltonian, with 
the cosine interaction at a point rather than along the line has appeared in various physical 

Note that the perturbation changes only a single Fourier 
Therefore in the limit L — > oo the perturbation has no effect, 



contexts, see for example Ref. |32| 
mode of the chiral fermion. 
which is to be expected since a term in the Hamiltonian consisting of a product of an odd 
number of fermions cannot produce any off-diagonal long range order. 



VII. SEMICLASSICAL CONSIDERATIONS 

In the case of the standard sine-Gordon theory, the classical equation of motion has 
massive solitonic solutions. A semiclassical analysis about these configurations, valid for 
small j3 2 , gives information about the spectrum of the quantum theory in the massive phase 



which exists for (5 < 4 [33]. In this section we demonstrate that the classical field theory 
of the xSG model possesses analogous solitary wave (kink) solutions, and we address the 
question as to whether or not a similar semiclassical expansion is useful. 

Consider the model (in the limit L — > oo) as a classical field theory. The real-time 
Lagrangian density is 

£ x sg = --rd x (/)(d t (j) + v n d x 4>) + k(cos(/30) - 1), (126) 

where the classical coupling constant k is related to the dimensionless coupling used in 
Section |T| by k = —A/ (27ra) 2 . We have included a constant term in the potential energy so 
that the minima at <ft = 2Tcr//3, r 6 Z have zero energy. The classical equation of motion is 

d x {d t + vM<f> - 2™/3sin(/30) = 0. (127) 

The conjugate momentum of the field <f>(x) is U(x) = —d x <p(x)/4:7r, and the Hamiltonian of 
the theory is 



n 



XSG 



J dx [Anv n (U(x)f - K(cos0(j)) - 1)] . (128) 



The equation of motion is reproduced in the Hamiltonian formalism by defining the funda- 
mental Poisson bracket of the chiral Bose field to be \EM 
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{d x (j)(x), 4>{x')}p B = 2tc5(x - x 



(129) 



This Poisson bracket is non-canonical, which is to be expected since the system is con- 
strained, i.e., Il(x) can be expressed in terms of <p(x). This is in contrast to the standard 
SG theory, in which the field and its conjugate momentum are independent dynamical vari- 
ables. If we define the rescaled field = f3(p to normalize the period of the potential energy 
term, we see from Eq. (|129| ) that the Poisson bracket of d x <p(x) and 4>{x') is proportional to 
p 2 . Therefore we would expect a semiclassical analysis to be relevant when p 2 is small, just 
as in the non-chiral SG model. However, recall that the quantum xSG theory at a finite 
length L is only sensible for p 1 e Z, and hence in the region of p 1 where we expect the 
semiclassical expansion to be most valid, the quantum theory is problematic. 

If we consider only static field configurations, the equation of motion is identical to the 
one obtained from the non-chiral SG theory. Therefore, Eq. 



]127| ) has time-independent 
kink solutions which are identical to those of the standard SG theory. From these stationary 
configurations we can find time-dependent solutions by applying a boost. The result is: 



<f>±,v a {t,x) 
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(x - v s t) 



(130) 



where v s E [— oo,v n ] is the velocity of the solitary wave. The topological charge is 

N = ^- [ dxd x cj)(x), 

Z7T J 



(131) 



and the solution (j)± tVs has N = ±1. These solitary waves are not chiral, v s can be positive 
or negative, but there still is left-right asymmetry. The solitary waves that move to the 
right (the chirality at k — 0) have a maximum velocity v n and are contracted relative to 
the v s = solutions, while the solitary waves moving to the left have no maximum velocity 
and are stretched relative to the static solutions. The existence of solitary wave solutions 
with arbitrarily large velocities is in contrast to the standard SG theory where the solitons 
have a maximum velocity because of the Lorentz invariance of the theory. Also, in the usual 
SG theory the solitary wave solutions are solitons. Whether or not this is true in the %SG 
theory is an open question. 

Since C x sg does not depend explicitly on t or x there is a conserved energy, £, and a 
conserved momentum, V, given by 



dx 



4:7V 



{dAf - /<cos(/30) - 1) 



n<P]= I dx— {d^f 



(132) 
(133) 



Using Eq. fll30| ) in Eqns. (|132|) and ( |133| ), the energy and momentum of the solitary wave 
solutions are 



E(v s ) = E[(t>±,v a . 
P(v s ) = V[^ Va 



VnC 



c 

V V n - Vs 



(134) 
(135) 
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where ( = \J8k/7t(3 2 . Note that P(v s ) > even for solutions that propagate to the left (v a < 
0). The dispersion relation is E(P) = v n P + ( 2 /P, which has a minimum at P = C/s/^n 
for which E(P Q ) = 2^(. 

We have found that there are topologically non-trivial solitary wave solutions of the 
classical equation of motion which exhibit a dispersion relation with a finite energy gap of 
order y/R. However, the quantum spectrum of the x^G theory is gap less, at least at the 
points p 2 = 2 and p 2 = 4. It is unclear how a semiclassical expansion about massive classical 
solutions could possibly describe the massless quantum theory. A second observation is that 
while the solitary wave solutions in Eq. ( |130| ) conserve the topological charge, the dynamics 
of the xSG theory does not. From Eqns. ( |128| ), ( |129|) , and (|131|) we see that 

{N, H xSG }pb = / dx sin P<f>(x). (136) 

The functional on the r.h.s. of Eq. ( |136| ) is in general non-zero, but it does vanish for <j)± tVa 
since they are odd functions of (x — v s t). If at some instant an arbitrarily small perturbation 
which destroys this symmetry is added to the solitary wave solution, it follows from Eq. (|136|) 
that a naive time evolution would result in a non-conserved winding number. However, from 
Eqns. ( |131| ) and ( |132| ), we see that if the topological charge of the system departs from an 
integer value, the energy would be infinite. Thus, if we impose boundary conditions on the 
field and its spatial derivative such that the energy is finite, specifically, 

lim d x 6(x) = 0, and lim 6(x) G f ^Z- ) Z, (137) 

x— »±oo x— >±oo \ J 

then the winding number will be restricted to the integers, and hence it will be conserved 
if we assume continuity of <p(t,x). In the quantum %SG theory the topological charge is 
not conserved, as can be seen from the non-vanishing of the commutator corresponding 
to the Poisson bracket in Eq. ( |136| ). We therefore believe a semiclassical expansion about 
the solitary wave solutions is unlikely to be useful for understanding the spectrum of the 
quantum theory. 



VIII. SUMMARY 

We have considered the edge theory of double-layer quantum Hall systems with a single 
mode per edge, upon inclusion of interlayer single-electron tunneling at the edge. The theory 
can be separated into a free chiral boson Hamiltonian for the charged mode and a chiral 
sine-Gordon Hamiltonian for the neutral mode. We have argued that the RG flow of the 
xSG theory differs significantly from that of the standard, non-chiral sine-Gordon theory, in 
particular (3 is not renormalized and there is a line of fixed points at p 2 = 4. In addition, 
the xSG theory likely only makes sense for integer values of p 2 . 

These values in turn break up into an apparently unphysical fermionic sequence, and a 
bosonic sequence that is realized in quantum Hall systems. In the bosonic sequence there 
are two points at which tunneling is not irrelevant, p 2 = 2 and p 2 = 4, and each corresponds 
to an infinite set of double-layer quantum Hall edge theories. For these two cases we have 
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obtained exact solutions for the partition function and some correlation functions. At p 1 = 2 
we find that tunneling produces spatial oscillations in the correlation functions while at 
p 1 = 4 tunneling splits the velocities of the two Majorana components of the chiral Dirac 
fermion corresponding to the neutral mode. 

The exact solution of the p 1 = 4 case involves an unfamiliar Bose- Fermi identity, ipd x ip ~ 
e 2 ^, which is put on firm ground in the appendices. This is the first in a sequence of identities 
of the form 

k-l 

e m*) „ TJ d™tp(x), k e Z + . (138) 

m=0 

Extensions of this work to include the effects of disorder, asymmetric velocities v ^ v', 
and the case of more than two layers are currently being pursued and we expect to discuss 
these, along with a proposal to detect the trifurcation of the electron at the 331 edge, in a 
future publication. 
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literature we have chosen to retain our proofs in the interests of completeness. 

APPENDIX A: KLEIN FACTORS 

In this appendix we demonstrate that it is possible to modify our definitions of fermionic 
vertex operators in a way that gives the proper anticommutation relations between different 
flavors without altering the form of the Hamiltonian. The first case to consider is the physical 
electron operators defined in Eq. (f25|), which commute for electrons in different layers. We 
modify this definition to 

= j^e™ s ^e lu ^\ (Al) 

where S is a matrix used to enforce the correct anticommutation relations. We now have, 
using Eqns. ffi) and fllPp, 

* i (x)* i (ar') = ^^x'^^xy^mix'-x^-s^K^+s^Kui) for x ^ x \ (A2) 
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If K12 = n is odd we can set S = and have the correct anticommutation relations between 
the electron operators in different layers. If n is even then we require 

n(S n - S 22 ) + m'Sn - mS 2 i E Z odd . (A3) 

From the dimension of the operator in Eq. (|A1|) it is clear that for m 7^ 1 it is not possible 
for *&i(x) and ^\(x') to have a canonical anticommutator; a characteristic consequence of a 
low-energy projection. It can be shown that {: *&i(x) :, : ^(x 1 ) :} oc 5ij5 <K ' m ~ l \x — x') in the 
limit a — > 0, L — > 00. 



Hl = — [ L/2 dx ( :e i ^2j(K r 2~K ]1 )-(S l3 ~S 23 W J ) e iu 1 (x) e -iu 2 (x) . + h( o ( A4 ) 
L m J-L/2 



Using Eq. ( |A1|) in the definition of the tunneling Hamiltonian (29|) gives: 

An r L / 2 

:/2 

The topological-charge dependent coefficient is just ±1 provided: 

5*12 — 5*22, 5*11 — 5*21 E Z 

even • 

(A5) 

For any even n and odd m, m' there are infinitely many choices for S that simultaneously 
satisfy Eq. (|A3|) and Eq. (A5). We make the particularly simple choice 

S=(\ °) for neZ even , (A6) 

and absorb the overall sign into the definition of the tunneling amplitude A — ► (— l) n+1 A, 
making Eq. ( |A4| ) identical to Eq. 

The second place where Klein factors are needed is in the fermionizations of the 
Hamiltonian in Section |V[ For the case p 1 = 2, Eq. (R^) should be modified to read 

tfjiix) = —, =e i * Si * N i e**® . (A7) 
V27ra 



The matrix S is again given by Eq. (|A6|) and ensures the proper anticommutation relations 
between the different fermion species without modifying Eq. just as in Eq. (|A4|) . For 
the multi-flavor fermionization at p 1 = 4, Eq. (|105|) should be replaced by 

if>i(x) = -^e^^e^W, (A8) 
\/2na 

where the matrix P is chosen so that ipi and ipj anticommute for i ^ j. One can readily 
show: 

^(x)^i(x') = Vi(^ / )^(^)e i7r(Pii ~ Py+sgIl(a: ~ a: ' ),5 " ) for x^x', (A9) 
and therefore we demand 

Pij-PjiEZow for i^j. (A10) 
Using Eq. ( |A5| ) in Eq. ( |103| ) we find the fermionized Hamiltonian including Klein factors is 
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n 



XSG 



L/2 
-L/2 



dx : 



iv'J)\d x ipi + A(e 



iirp inOiPijN. 



^^lWa^i + h.c. 



(All) 



where O t = (-1) 4 and p = P 21 - P 22 + P32 ~ P31 + Pa ~ P42 + P43 ~ Pa- From Eqns. flAToD 
and (|A11 ) we see that one choice of P that gives the correct anticommutation relations and 
which trivializes the explicit topological charge dependence of the Hamiltonian is 



P 







1 1 ix 

1 1 1 

1 

Vo 1/ 



(A12) 



With this choice of P, we find Eq. ([All] ) is identical to Eq. (|106| ). 



APPENDIX B: EQUIVALENCE OF CORRELATION FUNCTIONS 

In this appendix we prove the validity of the identity 



il)(x)d x i/j(x) 



2vra 2 



2i4,{x) 



(Bl) 



by comparing the 2iV-point correlation functions of these operators and their Hermitian 
conjugates. We work at zero temperature for simplicity Consider first the correlation 
function in the free chiral boson theory: 



N 



B N (z h . . . , z N ; 4, • • • , z' N ) = ( Yl 



\i=l 



2na 2 



n 

i=l 



\2na 2 J 



1 Ili<i<j<N( z i — z j) 4 { z 'i ~ z 'j^ 4 



{2k] 



2N 



T\l<i,j<N( z i ~ z j) 4 



(B2) 



From Eq. (|B1|), the corresponding correlation function in the free chiral fermion theory is: 

IN N \ 

T N {zi, . . . , z N ; z[, . . . , z' N ) = m if)(ri, Xi)d Xi i){Ti, x t ) Y[ {~^{ T i x 'i) d x'^H T i X 'S) ) • ( B3 ) 



We adopt the complex notation ip(r,x) = tp(z),d x ip(r,x) = —id z ip(z). To evaluate this 
2A^-point function we write it as the limit of a 4A^-point function and pull the derivatives 
outside the correlator: 



FN&Zi}, i z i}) = . , li, m . , lirn 



lim 



lim 



W 1 ,W 2 ^z 1 w' v w' 2 -*z[ W2N-1>V>2N-*ZN W 2N-1' W 2N^ Z 'n 



(B4) 



j=l \i=l i=l 

From Eqns. ( PD| ) and Q5T]) we find that the 4A^-point correlation function can be written: 
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II V» W II ^K) = TnZvw rr r.„. _ „,n ' ( ) 



\i=i j=i 



(27T)^ ni< M <27vK-^) 



We first organize the products appearing in Eq. (|B5| ) by separating coordinates, Wj, with 
even and odd indices. The denominator can be expressed as 



2N 2N N N 



n ~ w 'j) = n n [(^-i - w' 2j _ 1 )(w 2i - u^-ix^-i - w 2j ){w 2i - u>y] > ( Be ) 

i=l j'=l i=l j'=l 

and the numerator can be rewritten using 

2AT-1 2N N-l N r 

n n ~ w i) = n n (^a-i - ^-lX^ - 

i=l j'=i-f 1 i=l j=i+l 

1 ^ 

x(w 2 i-i - w 2 j)(w 2 i - w 2 j) Y[(w 2k ^i - w 2 k), (B7) 

J fe=i 

and an identical expression with u>j replaced by w\. 

When the N 2 derivatives in Eq. ( [B4| ) act on the expression on the r.h.s. of Eq. (|B5|) , a 
large number of terms are generated. However, from Eq. ( |B7|) we see that the numerator in 
Eq. ( PS] ) contains a factor (w 2 k-i — w 2 k) for each k = 1, . . . N. In the limit w 2 k-i, w 2 k — > Zk, 
this factor vanishes, and hence the only term that will be non-zero in the multiple limit 
given in Eq. (|B~|) is the term in which the operator d W2k acts on {w^k-i — w 2k ) for all 
k = 1, . . . N. An analogous argument holds for the primed coordinates and therefore we find 
using Eqns. (jBSJ), and (g|-([B7p: 



•^jvd^il, {z'A) = lim lim ... lim lim 

w u w 2 -*zi w' 1 ,w' 2 ^z' 1 u>2N-\,w 2N ->zn w' 2N _ 1 ,w' 2N ^z' n 
1 IIl<i<j<iV [(W 2 j-1 ~ W 2 j-l)(w 2 i - W 2 j-l)(w 2 i-l - W 2j )(w 2i - W 2j )} X [Wj -» w'jl 
(2tt) 2N Ul<i,j<N(w2i-l - W2j_l)(w 2 i - W , 2j _ 1 )(w 2 i-1 ~ W' 2j )(w 2i - w' 2j ) 

< i<3 <N(z t - Zj n4-z- _ Bif{{Zi}M})i (B8) 



1 Hl<i<j<N( z i ~ z j) 4 ( z i ~ z '^ ' 

WF rii< M -<v(^ - z \ 



completing the proof. 



APPENDIX C: EQUIVALENCE OF COMMUTATION RELATIONS 

Here we calculate some additional commutation relations in support of the identifica- 
tion ([BSD, which can be written in terms of a normal-ordered exponential [using Eq. (|T0|)] 

as 

27T? 

— :e**: = ^. (CI) 

Consider the operators 
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V K \ n 

C^\x) = - : (e 2i * ± e" 2 ^) : = ]T e ik ^ 



(C2) 
(C3) 



where hi = 2n/L and the rightmost sides of the expressions define the normalization of the 
associated harmonics. One can check that 



d(±) a(°) 

P(±) A (1) 



-2C 



(T) 
n+mi 



in — m) C 



"r;(±) p(±) 
C (+) C H 

For example, to derive Eq. ( 



(±) 

n+mi 
(1) 



(n - m) A^| m + — n (n 2 - 1) 5 n+m , 
(2) n 2 -nm + m 2 -l (Q) 



we write each operator clS cL Fourier transformation, 

i-L/2 



(C4) 
(C5) 

(C6) 
(C7) 



CW = L" 1 / dxe- i * wX C (±) i 

i-L/2 



;r 



and evaluate the products of normal ordered vertex operators with the help of decomposi- 
tion ([[!]) explicitly at a given cutoff a > 0, 



. Uxx<j>{x) . . e ia 2 <j>(y) . _ . iai<j>(x)+ia2<l>(y) 



-iriaia'2 (y—x)/L 



I — e 2iri(y-x+ia)/L 



(C8) 



In the limit a — > 0, the commutator of these two vertex operators vanishes at x ^ y for even 
integer values of the product «i«2, but it has a pole of order (aia^l at x = y for negative 
values of the product, a% a% < 0. A direct calculation with the help of the substitution 
z = exp(27r % x I L) shows that the only contribution to the commutator ( |C6| ) is given by the 
corresponding residue, 



p(+) p(+) 



1 



E 



a ~ ±2 \ z \=l v ' v y 



a 2 /2-m-l . gia (<p z —<Pw) . 



(C9) 



24 



E 

a=±2 



1-n 



5 3 



27T2 9w 3 



UI = 2 

d 2±riL pS^l + {m - 1) * : ( W : - ^ : ( W : 



1*1=1 



12 v r 7 2 

and a simple integration by parts gives the result ( |C6| ) . In Eq. ( |C10| ) we used [cf. Eq. (|1 

oo -i 

Vz = -NR ilnz-x/R + J2^p e ~ kra/2 { z ' r b t + z ' r b r) ■ 



(CIO) 
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Readers familiar with conformal field theory techniques will notice that Eq. ( |C9|) could be 
obtained immediately within the formalism of canonical quantization on a cylinder [pER . We 
also notice that Eq. ( |CJ5| ) follows from a more general relationship for the vertex operator 
V a (x) =:e iQ ^:, 



V a A W 



(Cll) 



which can be derived in a variety of standard ways. 

On the other hand, using Eq. ( |CT| ) and the canonical bosonization prescription, we can 
define the fermionic counterparts of the same operators, 



A (p_1) (x) = — :{:c c:)":. 

p \ ) 

C ^ (x) = ^(^ W T ^ V f ) = 5 E e * ^ E s {cr-c. ± 



where the canonical Fermi operators c r , {c r , cj} = <5 rs are defined by the expansion of the 
ant iper iodic field 

V(x) = L- 1 / 2 Yl c r e 2mrx/L . 

r6Z+I 

Comparing their commutation relations, it is a straightforward (but lengthy) exercise to 
express the operators 

a(0) — ft (0,0) A (1) _ 1 f -6 (1,0) , -d (0,1) \ A ( 2 ) — f> (1,1) i + 1 f> (0,0) 



in terms of operators 



K+n 

Bfr*= lim £ (r-n)V 6 :cj_ n c r :, 

if— >oo 



bilinear in fermions. The commutation relations then read 



fl(±) iM 
p(±) 



c(±) r(±) 
c(+) c(-) 



-2C 



(T) 
n+mi 



(±) 



(n — m) C 
(n - m) AJ2 m + (n 2 - 1) 5 n+i 
2 pd.D | (n + m) 2 -nm ~ (0 , 0) 

~(2) n 2 -nm + m 2 -l - (Q) 



6 



they are in precise correspondence with Eqns. ( |C4| )-(P7|). 
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APPENDIX D: ALTERNATIVE SOLUTION OF THE xSG THEORY AT (3 2 = 4. 



With the help of the commutation relations derived in Appendix |C], we can diagonalize 
the xSG Hamiltonian at p 2 = 4, Eq. flS5|), directly within the operator formalism. Indeed, 
in terms of normal-ordered quantities the Hamiltonian reads 



n 



(3 2 =4 
X SG 



L/2 

dx 

-L/2 

2tt r L / 2 



^:(^0) 2 :+^:cos(20): 



L 2 J-L/2 



dx 



A 



v n A (1 \x) + -C {+ \x) 

71 



2tt 

T 



7T 



(+) 



(Dl) 
(D2) 

where we used the normal ordering formula (0) and the definitions ( |C2j ) , ( |C3| ) ; the last 
line is written as a sum of two commuting operators. The commutation relations (|C4] )- (|C7|) 
imply that the operators 

generate two independent c = 1/2 Virasoro algebras, 



j a Tb 

1J ni m 



5 ab (n - m) L b n+m + ±n (n 2 - 1) <W,o S ab 



and the Hamiltonian ( [D2|) becomes 

^XSG = ( 2n / L ) ( V + L + V - L ) 7 V ± = v n± A/VT. 



Now the spectrum can be obtained directly in terms of irreducible representations of the 
two c = 1/2 Virasoro algebras. To this end it is convenient to write each set of operators 
a = ±, in terms of an independent Majorana fermion. In particular, 



2tt 1 f L / 2 

~T L o = 7T / dx : Xad x Xa ■ 

L 2l J-L/2 



which takes us back to the Hamiltonian ( |70| ) obtained by a more conventional fermionization 
procedure. 
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